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Abstract — In this paper, theoretical limits on time delay 
estimation are studied for ultra-wideband (UWB) cognitive radio 
systems. For a generic UWB spectrum with dispersed bands, 
the Cramer-Rao lower bound (CRLB) is derived for unknown 
channel coefficients and carrier-frequency offsets (CFOs). Then, 
the effects of unknown channel coefficients and CFOs are 
investigated for linearly and non-linearly modulated training 
signals by obtaining specific CRLB expressions. It is shown that 
for linear modulations with a constant envelope, the effects of 
the unknown parameters can be mitigated. Finally, numerical 
results, which support the theoretical analysis, are presented. 

Index Terms — Ultra-wideband (UWB) cognitive radio, time 
delay estimation, Cramer-Rao lower bound (CRLB), carrier 
frequency offset (CFO), dispersed spectrum utilization. 

I. Introduction 

Cognitive radio systems provide intelligent wireless com- 
munications by means of adaptation, feedback, learning, and 
spectrum sensing. The main characteristics of cognitive radio 
are related to spectrum utilization in a dispersed manner; i.e., 
a cognitive radio system can utilize a number of adjacent or 
non-adjacent frequency bands at the same time JTJ , 0. In an 
ultra-wideband (UWB) cognitive radio system, the available 
spectrum has ultra- wide bandwidth 0, and the spectrum can 
be used in a dispersed manner J4), Q. 

Cognitive radio systems with location awareness provide 
location information to the users, which facilitate location- 
based services and applications, and also help network op- 
timization J6), 0. The part of the cognitive radio system 
that provides location sensing is called a cognitive positioning 
system (CPS) J8). CPSs can use Cramer-Rao lower bound 
(CRLB) information at the transmitter side to optimize the 
transmission parameters for achieving specific accuracy re- 
quirements. In 0, the CRLB for time delay estimation in 
multipath channels is derived, where the path delays and 
coefficients are assumed to be unknown. An extension of 
J8) is provided in Q)), which obtains the frequency-domain 
CRLB in dynamic spectrum access systems for unknown path 
delays, distance-dependent coefficients, phases and frequency- 
dependent coefficients. Both J8) and J9) consider the available 
bandwidth in the form of a single piece (i.e., the whole band- 
width). However, the available bandwidths in UWB cognitive 
radio systems are commonly dispersed ifTOl and each user 
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can transmit and receive over multiple dispersed bands. This 
study considers UWB cognitive radio systems that can employ 
various dispersed bands IfTTl . 

In this paper, CRLBs for time delay estimation are studied 
for UWB cognitive radio systems with dispersed frequency 
bands, and the effects of unknown channel coefficients and 
carrier frequency offsets (CFOs) are quantified. After ob- 
taining a generic expression, specific CRLB formulas are 
obtained for various modulation schemes. Then, it is shown 
that the same lower bound can be achieved for the cases of 
known and unknown CFOs for linear modulation schemes. 
In addition, it is proven that the effects of unknown channel 
coefficients can be mitigated for linear modulations with a 
constant envelope. Finally, numerical results are presented and 
concluding remarks are made. 

II. Signal Model 

A UWB cognitive radio system that occupies K different 
frequency bands, which can be adjacent or non-adjacent, is 
considered. The transmitter sends a signal occupying all the 
K bands, and the receiver tries to estimate the time delay 
of the incoming signal. Such systems can be implemented as 
orthogonal frequency division multiplexing (OFDM) systems, 
where the sub-carriers corresponding to the unused bands 
are considered to have zero coefficients lfT2l . |[T3"1 . The main 
disadvantage of this approach is that when the available spec- 
trum is very dispersed, processing of very large bandwidths 
is required, which can increase the receiver complexity IfTTl . 
lfT4l . Another approach to implement cognitive radio systems 
with dispersed spectrum utilization is to process the signal in 
multiple branches, as shown in Fig. fTl IfTTl . In this case, each 
branch processes one available band, which is significantly 
narrower than the total bandwidth. Unlike the first scenario, the 
system uses multiple down-conversion units in Fig. [T] hence, 
different carrier frequency offsets (CFOs) can be observed 
at the signals at different branches. Considering the receiver 
structure in Fig. \T\ the aim is to determine the accuracy of time 
delay estimation in the presence of unknown channel effects 
and CFOs. 

Since each of the K frequency bands will have narrow 
bandwidths compared to the total system bandwidth, it is 
assumed that the baseband representation of the received signal 
in the ith branch can be expressed as 

n{t) = a^-^s^t - r) + m(t) , (1) 
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Fig. 1 

Block diagram of a cognitive radio receiver with multiple 
branches. 



for i = l,...,K, where 0*4 = aiC?'^ i and u>i represent, 
respectively, the channel coefficient and the CFO for the signal 
in the ith branch, Si(i) is the baseband representation of the 
transmitted signal corresponding to the ith band, r is the time 
delay, and rii(t) is complex Gaussian noise with independent 
and white components, each having spectral density of. 



III. Generic Limits 



In the presence of unknown channel and CFOs, the vector 
of unknown signal parameters can be expressed as 6 = 
[t ai ■ ■ ■ a,K <pi ■ ■ ■ 4>k 0J1 ■ ■ ■ ujk]- Then, for an observation 
interval of [0,T], the log-likelihood function for can be 
obtained from dTJ as 
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where c is a constant that is independent of 6. 

In order to calculate the CRLB on time delay estimation, 
the Fisher information matrix (FIM) should be calculated first 
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diagonal matrix with its ith diagonal being equal to Xi, Ei = 
J \si(t — T)\ 2 dt is the signal energjQ, Ei is the energy of the 
first derivative of s,(i); i.e., Ei = |s^(i — t)| 2 dt, and Ef-, 
E], Fi, F and G, are given, respectively, b>Q 

rT 
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The CRLB for unbiased time delay estimators can be 
obtained from the element in the first row and first column 
of the inverse of the FIM in (0, i.e., [I -1 ] . Based on the 
formulas for block matrix inversion, it can be calculated as 

Ell 
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where 
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In order to investigate the effects of unknown CFOs, the 

'Although Ei is a function of r in general, it is not shown explicitly for 
convenience. The same convention is used in U2M16I as well. 
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CRLB for time delay estimation can be calculated in the 
presence of known CFOs. In that case, similar derivations 
can be performed for the unknown parameter vector 6 = 
[r fli ■ ■ ■ clk 0i • • ■ <t>K\, and the CRLB can be obtained as 

CRLB 2 = — j^- : , (19) 
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where 
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Note that CRLB 2 in dT9jl is smaller than or equal to CRLBi 
in dTD i in general, as more unknown parameters exist in the 
latter case. 

In the extreme case in which both the channel coefficients 
and the CFOs are known, the unknown parameter vector 
reduces to r, and the CRLB expression reduces to 

CRLB 3 = 1 . (21) 

Ei=l H E i 

IV. Special Cases 

In this section, the baseband signal s 2 (i) in ([TJ is modeled 
to consist of a sequence of modulated pulses; i.e., 



i (t)=J2d i ,i(t)p i {t-lT i ) 



(22) 



for i = 1,...,K, where di t i(t) denotes the complex data^ 
for the Zth symbol of signal i, and Pi(t) represents a pulse 
with duration Tj, which is non-zero only for t € [0,2}]. In 
addition, it is assumed for the simplicity of the expressions 
that the observation interval T can be expressed as T — NiTi 
for an integer Ni for i = 1, . . . , K. 

First, consider a scenario in which the signal Sj(i) is a 
linearly modulated signal expressed as Si(t) = Yli di,iPi(t — 
ITi). In this case, it can be shown from ( T22l that Sj(t)s|(t) 
is a real quantity; hence, CRLBi in ( fT7] i and CRLB 2 in ( fT9l 
become ifTTl 

CRLBi = CRLB 2 = — — - ■ (23) 
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The result in ( f23l l implies that for common modulation types 
such as pulse amplitude modulation (PAM), phase shift keying 
(PSK) and quadrature amplitude modulation (QAM) |16|, the 
CRLB of time delay estimation for the case of unknown CFOs 
is the same as that for the case of known CFOs. 

Next, consider linearly modulated pulses with a constant 
envelope; i.e., Si(t) = J2i <k,lPi(t - ITi), with \d it i\ = \di\ VZ. 
In addition, if Pi(t) satisfies p,(0) = Pi(Ti) for i = 1, . . . ,K, 
the CRLBs in ( fTTI i and ( fT9l can be calculated as 



CRLBi = CRLB 2 = CRLB 3 = 



^ELsnr^ 



, (24) 



where SNR 4 = NAdA 21 - 



A 2 e Pi 



with E Pi = J^pKtfdt, and 



From ( [24-b . it is observed that for linear modulations with a 
constant envelope, such as PSK, the CRLB on time delay esti- 
mation is the same for known or unknown CFOs and/or chan- 
nel coefficients when the pulses satisfy p(0) = p(T), which is 
commonly the case in practice. Therefore, in such cases, it is 
possible to mitigate the effects of the unknown parameters on 
time delay estimation. More specifically, maximum likelihood 
(ML) estimators of time delay can asymptotically achieve the 
same CRLB for unknown CFOs and channel coefficients case 
as the ones in presence of CFO and/or channel coefficient 
information lfl5l . 

V. Numerical Results and Conclusions 

In this section, numerical studies are performed in order to 
evaluate the CRLB expressions in the previous sections. For 
the cognitive system, it is assumed that all the K bands in 
the system have the same bandwidth and the same pulse is 
employed for all of them; i.e., pi(t) = p(t) for i = 1, . . . , K 
(c.f. d22l). For the pulse shape, the following Gaussian doublet 
is employed 
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Pi is the effective bandwidth of pi(t) IfTTl . 



where A and £ are parameters that are used to adjust the pulse 
energy and the pulse width, respectively. In the following, 
£ = 8 ns is employed, for which the pulse width becomes 
approximately 20 ns, and A is selected to normalize the pulse 
energy. In addition, the same noise spectral density is assumed 
for all the K branches; i.e., <7j = a for i = 1, . . . , K. Finally, 
the system SNR is defined as the sum of the SNRs in the 
various branches. 

In Fig. [2] the CRLB expressions in (fTT), (QUi and (ETT) . 
which are labeled as CRLBi, CRLB2 and CRLB3, respec- 
tively, are plotted versus SNR for three different modulation 
types, namely 16FSKQ, 16QAM, and 16PSK. In all cases, the 
same modulation sequence is employed at different branches; 
i.e., dij(t) = di(t) for i = 1,...,K, and that the channel 
amplitudes are normalized to unity; i.e., |a»| = 1 for i = 
1, . . . , K. In addition, the modulation sequence for each mod- 
ulation type is scaled appropriately so as to equate the energy 
of the sequences in different scenarios. Also, there are K = 10 
branches in the system, and N — 2 symbols are received 
at each branch. It is observed that for all the modulations, 
CRLB 3 < CRLB 2 < CRLBi is satisfied, since CRLBi 
corresponds to the case of unknown delay, channel coefficients 
and CFOs, CRLB 2 corresponds to the case of unknown delay 
and channel coefficients, and CRLB3 corresponds to the case 
of unknown delay only. In other words, for the cases with 
fewer unknown parameters, lower CRLBs are observed. For 
the 16FSK modulation, all three bounds are distinct. Note that 
as 16FSK is a non-linear modulation, the result in ( |23l l does 
not apply in this case. For the 16QAM case, CRLBi = CRLB 2 
as expected from d23l l. However, CRLB 3 is lower than both, 
which is can happen since 16QAM is not a constant envelope 
modulation (c.f. (1241). Finally, for the 16PSK case, all the 
three bounds are equal, which verifies (f24l . 

In Fig. [3] N = 16 symbols are used and the rest of the 
system parameters are kept the same as in the previous case. 



3 Since data-aided time delay estimation is considered, known data symbols 
are assumed. 



4 The amount of frequency shift is selected as 23.75 MHz for the FSK 
modulation. 
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VCRLB VERSUS SNR FOR K = 10 AND N = 2. 



Fig. 3 



V / CRLB VERSUS SNR FOR K = 10 AND N = 16. 



It is observed that the results are similar to those in Fig. [2] 
except that the gap between the CRLBs is decreased, and lower 
CRLBs are achieved when more symbols are employed. In 
addition, the reduction of the gaps between CRLBi, CRLB2 
and CRLB3 implies that using a larger number of observation 
symbols can mitigate the effects of unknown CFOs and/or 
channel coefficients. 

Finally, the CRLBs for different numbers of branches, which 
represent the numbers of available dispersed bands in the 
spectrum, are investigated for various SNRs. The parameters 
are the same as in the previous scenario, except that only 
16PSK is considered here for simplicity, and the CRLB is 
evaluated for K = 10, 15, 20, 25 branches and for SNR = 
5, 10, 15 dB. It is observed from Fig. that the CRLBs 
decrease as the SNR increases but they are invariant to the 
changes in the number of branches. This is because the SNR 
is defined as the sum of the SNRs in the different branches; 
hence, the SNR per branch is reduced as more branches are 
employed. 
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